Top Quark Pairs at High Invariant Mass — 
A Model-Independent Discriminator of New Physics at the LHC 
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We study top quark pair production to probe new physics at the LHC. We propose reconstruction 
methods for ti semileptonic events and use them to reconstruct the ti invariant mass. The angular 
distribution of top quarks in their cm. frame can determine the spin and production subprocess for 
each new physics resonance. Forward-backward asymmetry and CP-odd variables can be constructed 
to further delineate the nature of new physics. We parametrize the new resonances with a few generic 
parameters and show high invariant mass top pair production may provide an early indicator for 
new physics beyond the Standard Model. 



In the next few years, high energy physics will experi- 
ence the excitement of major discoveries when the CERN 
Large Hadron Collider (LHC) opens up the unexplored 
TeV energy scale. Besides the long anticipated Higgs 
boson that is responsible for the mass generation in the 
highly successful Standard Model (SM), the quadratic 
sensitivity of the Higgs boson mass to radiative correc- 
tions indicate the existence of new physics associated 
with electroweak symmetry breaking naturally at the 
scale of order Airv. 1 Numerous extensions to the Stan- 
dard Model (SM) have been proposed to describe elec- 
troweak symmetry breaking. A sample of popular sce- 
narios include: the Minimal Supersymmetric Standard 
Model (MSSM) [K], new strong dynamics @, 0, H i], 
composite Higgs at the TeV scale |6j , Little Higgs theory 
fz[, and extra dimensions at the electroweak scale 
String-inspired extensions in the gauge sector associated 
with an extra U(l) symmetry [10j also lead to striking 
signatures. It is therefore highly expected that many new 
signatures will become manifest at TeV energy scales that 
can be probed at the LHC. 

The LHC will be a "top factory": About 80 million 
ti events will be generated by QCD production with an 
integrated luminosity of 100 fb _1 . Thus studying top- 
quark production can be fruitful. The fact that the top 
quark mass is at the electroweak scale (m t w v/y/2) 
suggests that top-quark production may be sensitive to 
new physics near the TeV scale. Generically, if the new 
physics contributes to ti production as an s-channel res- 
onance, we want to identify the signal as a bump on 
the smoothly falling ti invariant mass distribution. Once 
we can reconstruct the ti cm. frame, the integer spin 
(J = 0,1,2) of any resonances can be determined from 
the polar angular distribution of the top quark. An 
asymmetry of this distribution would probe the chiral 
structure of the couplings. It may be possible to explore 
the CP property of the couplings with the help of CP- 
odd kinematical variables constructed from the final state 



particle momenta. Moreover, the relative importance of 
gluon-gluon and quark-antiquark subprocesses can be in- 
ferred from the spin and angular distributions. It is thus 
of fundamental importance to effectively reconstruct the 

invariant mass via their decay products. 

We focus on the semileptonic decay mode, ti — > 
bj\j2 b£~u + c.c. where £ — e or /i. The purely hadronic 
decay mode of ti not only suffers from a much larger QCD 
background, but also loses the identification of t from i. 
For the purely leptonic mode, with a small branching 
fraction of about 4/81, one cannot reconstruct the ti in- 
variant mass with two missing neutrinos. The signal to 
search for is an isolated charged lepton plus missing en- 
ergy (ftr), 2 6-jets plus 2 light jets. The branching ratio 
of the semileptonic to the hadronic channel is 2/3. 

New Reconstruction Methods : 

A primary focus of our study is to reliably reconstruct 
ti kinematics at high invariant mass on an event-by-event 
basis. The challenge is to reconstruct the momentum of 
the missing neutrino. The transverse momentum of the 
neutrino is identified with the observed pr- The longitu- 
dinal momentum is subject to a two-fold ambiguity from 
solving the kinematic quadratic equation. 

Several top reconstruction methods have been used at 
the Tevatron 11 1. There, however, the top quarks are 



produced near threshold and the kinematics of the subse- 
quent decay products are very complicated. Since we are 
interested in new physics in the TeV region, demanding 
a high invariant mass for the ti events will tremendously 
simplify the kinematics, especially by distinguishing the 
b quark from b. Throughout our study, we use a 2 — > 6 
partonic level monte-carlo simulation that incorporates 
full spin correlations from production through decay [l2T ] . 
We made a Pythia simulation, including gluon radiation 
and hadronization, that confirmed our basic results. 

We impose a cluster transverse mass cut on the ti sys- 
tem 
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1 Here v ~ 246 GeV is the Higgs field vacuum expectation value. 



( Pb + Pi + Pn + Vn + pi) 2 + ft, + pr > 600 GeV. 

We adopt kinematical cuts from the ATLAS and CMS 
[l^ | top studies. We smear the hadronic energy ac- 
cording to a Gaussian error given by AEj/Ej — 
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0.5/ y/Ej /GeV © 0.03; and the lepton momentum by 
Ap T /p e T = 0M(p T /TeV) © O.O13/v / sin0, where 9 is the 
polar angle of the lepton with respect to the beam di- 
rection in the lab frame. We present two schemes to 
reconstruct semileptonic tt events and evaluate their ef- 
ficacy. 

(1). (Mw,rn t ) scheme: 

In this scheme, the key assumption is to take Mw and 
mt as inputs for their on-shcll production and decays. 



Step I: Demand 



The longitudinal momen- 



tum of the neutrino is formally expressed as 
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eT 



A PeL ± E t 



A 2 - Ap 2 



eT-i 



where A 



M w + 2p eT 



If A 2 - ApItE 2 > 0, 
the value of p V L that best yields the known top mass 
via m 2 ub — m 2 is selected. This ideal situation may not 
always hold when taking into account the detector reso- 
lutions. For cases with no real solutions, we then proceed 
to the next step. 

Step II: To better recover the correct kinematics, we in- 
stead first reconstruct the top quark directly by demand- 
ing mf ub = m 2 . The longitudinal momentum of the neu- 



trino is expressed as 



4 



PuL = A / p blL /2(E, 
x (p 2 blL A' 2 + (Ei 



pIil) ± 



PIil) 
4^H)) 1/2 , 



where A! 



M bl +2pur- fjr- The two-fold ambiguity 



is broken by choosing the value that best reconstructs 
= m 2 v . A plot of the top and W mass distributions 
is shown in Fig. 1(a). The solid histogram is from the 
procedure Step I, and the dashed histogram from Step 
II. With these two steps, there could still be some events 
that do not lead to a real solution, We thus discard them 
in our event collection. The discard rate is about 16%. 

(2). Small 9bw angle selection scheme: 2 
This scheme reconstructs the tt system without rely- 
ing on the top mass reconstruction, thus avoiding po- 
tentially large QCD corrections due to gluon radiations. 
Since we are interested in new physics in the TeV regime, 
the top quarks will be relativistic with a 7-factor of (1 
TeV)/2m t ~ 3. We thus expect that the decay products 
are fairly collimated along the top quark moving direc- 
tion. This is illustrated in Fig. 1(b) for the normalized 
opening-angle distribution between b and W + in the lab 
frame, where an increasing cluster transverse mass cut 



We have chosen to use O^yy in our reconstruction instead of, 
e.g., 6i v , because the b quark is on average much more energetic 
(highly boosted) than the lepton. 
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FIG. 1: (a) The W and top mass reconstructions from the 
{Mw, mt) scheme, with the procedure Step I (solid) and Step 
II (dashed), (b) Differential cross section versus cosOtw with 
no invariant mass cut (dotted), and with cuts of 600 GeV 
(dashed) and 1000 GeV (solid). 



has been imposed for the dotted, dashed and solid curves 
for M T > 0, 600, 1000 GeV, respectively. With the ini- 
tial requirement m 2 v — Mfy, the two- fold ambiguity in 
6bw is resolved by choosing the smaller reconstructed an- 
gle. This scheme should work better at higher tt invariant 
masses. 

Backgrounds to the tt Signal: 

The major backgrounds to our tt events include the 
processes W+ jets, Z+ jets, WW, WZ and ZZ. The 



ATLAS and CMS Technical Design Reports 14| detail 
studies of the selection efficiencies for these background 
processes in comparison to a reconstructed tt semilep- 
tonic signal. The ATLAS (CMS) group found for an in- 
tegrated luminosity of 10 fb _1 (1 fb _1 ) a signal to back- 
ground ratio of S/B = 65 (S/B = 26) 13]. Because 
of the expected high S/B ratio, our analysis is concen- 
trated solely on the tt events without including the small 
background contamination. Our analysis does not in- 
clude misidentification of faked leptons from jets in tt 
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FIG. 2: tt invariant mass distributions reconstructed from 
(a) the (Mw, mt) scheme, and (b) the small angle selection 
scheme. Both plots featured a 1 TeV resonance with a total 
width of 2% (solid), 5% (dashed), and 20% (dotted) of the 
resonance's mass. 

total hadronic decays. 

Although the t (t) is primarily identified by the charged 
lepton, £ + (£~), a concern is the matching of the b- 
jet associated with top quark decay. Both ATLAS and 
CMS studies [3] show a combination of kinematic fits, 
designed to properly reconstruct the W boson and the 
hadronically decaying top significantly reduces misiden- 
tification. Our cut on Mt helps significantly in this re- 
gard. 

Search for New Physics: 

We wish to explore the new physics searches in a 
model-independent manner for tt semi-leptonic decays. 
We consider tt production via 

gg -> <j) -> tt, qq-> V — >• tt, qq, gg — >• h -> tt, 

where <fi, V and h are the spin-0, spin-1, and spin-2 reso- 
nances. We characterize the effects on the invariant mass 
spectrum with three parameters: mass, total width, and 
the signal cross section normalization (w 2 ). The normal- 
ization u) = 1 defines our benchmark for the spin 0, 1 



and 2 resonances. They correspond to the SM-like Higgs 
boson, a Z 1 with electroweak coupling strength and left 
(L) or right (R) chiral couplings to SM fermions, and the 
Randall-Sundrum graviton h with the couplings scaled as 
A -1 for hqq, and (Ahx(Ml l /A))~ 1 for hgg, respectively. 3 
Numerically, we take A = 2 TeV. 

In Fig. 2 we show the reconstructed tt invariant mass 
distributions for the two reconstruction schemes. The 
SM tt total cross section is theoretically known beyond 
the leading order in QCD (l7f . We thus expect to have 
a good control of this distribution even at high invariant 
masses. As for new physics, we include the contribution 
of a 1 TeV vector resonance for illustration, for oj v — 
1, with total widths specified in the caption of Fig 2. 
We note that a very high invariant mass tail exists for 
the tt invariant mass reconstructed via the small angle 
selection. This comes from the mis-reconstructed events 
in the low invariant mass region. When a large enough 
transverse mass cut is applied for a given resonance, the 
tail will not obfuscate the resonance signal. For example, 
if searching for 2 TeV resonance, a 800 GeV minimum cut 
will eliminate the tail for the mass region of interest. 

We maximize the signal observability by isolating the 
resonance within an invariant mass window of ±100 GeV, 
±30 GeV and ±25 GeV for the scalar, vector and gravi- 
ton resonance, respectively. Given a resonance mass and 
total width, we can quantify how large to needs to be for 
a 5(7 discovery. With the number of events for a signal 
(S) and background (B), we require S/yB + S > 5. This 
translates to a bound uj 2 > (25 + 5\/25 + AB)/2S\ where 
Si is the benchmark signal rate for u> = 1. This is illus- 
trated by Fig. 3 versus the mass for a scalar, vector and 
graviton resonance for total widths of 20%, 5%, and 2% 
of its mass, respectively, for an integrated luminosity of 
10 fhr 1 . 

It is of critical importance to reconstruct the 
cm. frame of the resonant particle, where the funda- 
mental properties of the particle can be best studied. In 
Fig. 4, we show the top quark angular distribution, cos 8* , 
with 9* defined as the angle in the tt cm. frame be- 
tween the top-quark momentum and the incident quark 
momentum, with latter determined by the longitudinal 
boost direction of the cm. system. Although events in 
the forward and backward regions are suppressed due to 
the stringent kinematical cuts, we still see the impressive 
features of the d- function distributions [HI in Fig. 4(a): 
a flat distribution for a scalar resonance (dashed), d\ x 
distribution for the left/right chiral couplings of a vec- 
tor (dotted), and d\±± from qq (solid) and rf|±i from gg 
(dot-dashed) for a spin-2 resonance. To illustrate the sta- 
tistical sensitivity for observing the characteristic distri- 



3 More precisely, we use the Feynman rules given i n I l5l and in- 
clude the additional warp correction factors from |16|| . 
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FIG. 3: Normalization factor versus the resonance mass for 
the scalar (dashed) with a width-mass ratio of 20%, vector 
(dot-dashed) with 5%, and graviton (solid) 2%, respectively. 
The region above each curve represents values of ui that give 
5ct or greater statistical significance with 10 fb _1 integrated 
luminosity. 



butions, we show in Fig. 4(b) the expected SM tt events 
(solid) with lex statistical error bars in each bin for a 10 
fb _1 integrated luminosity, along with a 5<r signal of a 
chirally coupled vector summed with the tt background 
in the resonant region (dashed). Due to the large event 
sample, the statistical significance is evident in the cen- 
tral and forward region. The forward-backward asymme- 
try in cos 0* can thus be constructed to probe the chiral 
couplings of the particle to the top quark. With the 
identification of the charged leptons, one may even form 
kinematical triple products to test the CP properties of 
the couplings [19(. As for distinguishing the other spin 
states, it is a question of statistical significance. A nar- 
row graviton may be relatively easy to confirm, achiev- 
able with a few tens of fb _1 luminosity; a broad scalar 
may require the highest luminosity, more than 300 fb _1 . 

In summary, we investigated two ways to reconstruct 
semileptonic tt events at high tt invariant mass and 
showed the utility of each in discovering new physics in 
the form of integer-spin resonances. The angular distri- 
butions of the top in the reconstructed CM frame reveal 
the spin of the resonance, and relative contribution from 
the initial states qq or gg. The forward-backward asym- 
metry and CP-odd variables can be constructed to fur- 
ther differentiate models. Because SM top quark physics 
is well predicted, high invariant mass top pair production 
may provide an early indicator for new physics beyond 
the Standard Model at the LHC. 
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FIG. 4: Polar angular distributions for the top quark in the 
cm. frame, (a) Signal only by the (Mw, mt) scheme for a 
scalar (dashed), a vector (dots), and a graviton from qq (solid) 
or from gg (dot-dashed); (b) number of events for the SM 
tt background (solid) with la statistical error bars, and the 
background plus a vector resonance (dashed) . 
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